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ABSTRACT 



The Galactic population of open clusters provides an insight into star formation 
in the Galaxy. The open cluster catalogue by Dias et al. (2002b) is a rich source of 
data, including kinematic information. This large sample made it possible to carry 
out a systematic analysis of 481 open cluster orbits, using parameters based on orbit 
eccentricity and separation from the Galactic plane. These two parameters may be 
indicative of origin, and we find them to be correlated. We also find them to be 
correlated with metallicity, another parameter suggested elsewhere to be a marker 
for origin in that high values of any of these two parameters generally indicates a 
low metallicity ([Fe/H] Solar< —0.2 dex). The resulting analysis points to four open 
clusters in the catalogue being of extra-Galactic origin by impact of high velocity cloud 
on the disk: Berkeley21, 32, 99, and Melotte66, with a possible further four due to 
this origin (NGC2158, 2420, 7789, IC1311). A further three may be due to Galactic 
globular cluster impact on the disk i.e of internal Galactic origin (NGC6791, 1817, 
and 7044). 

Key words: galaxy: disk - open clusters: general. 



1 INTRODUCTION 

In general, Galactic open clusters are ensembles of stars 
with low central concentration and irregular shape. They 
constitute the outcome of a star forming region and are 
found across a range of evolving stages. Open clusters ar- 
guably represent the growth of the disk of the Galaxy via a 
smooth, self-regulated process of star formation, in contrast 
with globular clusters - formed during the early stages of 
the Milky Way - or with young massive clusters, found in 
starbursting galaxies. The fundamental difference between 
open and globular clusters appears to be the pressure of the 
environment, with open clusters formed in regions with low 
pressures (Elmegreen & Efremov 1997). The complex for- 
mation history of the Milky Way, characterized by a num- 
ber of events with enhanced star formation along with a 
smooth low-level constant formation (de la Fuente Marcos 
& de la Fuente Marcos 2004) should be imprinted in the 
"fossil record" represented by the open clusters, the glob- 
ular clusters and the bulk of the stellar population in the 
Galaxy. 

The total mass of an open cluster ranges from a few 
tens to a few thousands of solar masses (see for example 
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Durgapal & Pandey 2001), and have typical diameters of a 
few to a few tens of pc. Their stars form as dense clusters 
located within the parental giant molecular cloud (Pfalzner 
2009). Unlike the case of globular clusters, the population of 
open clusters is very strongly concentrated along the Galac- 
tic plane, which also makes them difficult to observe at large 
distances. As they evolve, open clusters undergo tidal evap- 
oration into the bulk of the stellar populations of the Galac- 
tic disk. This process is helped (Pfalzner 2009) if the clus- 
ter contains ~ 1,000 stars or more. With those numbers, a 
significant population of O stars is expected, whose stellar 
winds, together with supernovae from early-type stars, will 
expel the remaining gas, disturbing the clusters' dynami- 
cal equilibrium. Given this strong trend between cluster size 
and age, searches are often biased against old (i.e. more dis- 
persed) systems, harder to separate from the stellar back- 
ground. Hence, open clusters are usually found to contain 
bright, blue stars, indicating young ages. 

In addition to age, metallicity is a property that can re- 
veal the past history of stellar populations in clusters, disks 
and other components of the Galaxy. For example, Twarog 
et al. (1997) used metallicity to study the local structure 
and evolution of the disk, and Friel et al. (2002), used this 
property to determine the metallicity gradient of the disk. 
This type of study is also illustrated more recently by Fu 
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et al. (2009), and Magrini et al. (2009). The first group ex- 
amine the radial gradient of the metallicity of open clus- 
ters along the Galactic disk. From this they create a model 
whereby infalling gas grows the disk inside-out. Their model 
accounts for current and past metallicity radial trends. The 
second group find evidence for inflow that varies radially, 
being lower in the outer regions of the disk. This clearly 
establishes metallicity as a key parameter in deducing the 
past star formation history. 

There are about 2,000 known Galactic open clusters, 
but because they are difficult to observe, one should expect 
many more: Piskunov et al. (2006) estimate the total number 
of open clusters could be of order 10 5 . A large catalogue of 
open clusters was recently compiled from the literature by 
Dias et al. (2002b), henceforth referred to as DAML. The 
catalogue includes information about the kinematics of the 
clusters as well as estimates from the literature of their ages 
and metallicities. 

Several authors have used orbit calculations to study 
open clusters. Piskunov et al. (2006) considered the 
open cluster population in the all-sky catalogue ASCC- 
2.5(I/280A, Kharchenko, 2001) which is based on Tycho2, 
Hipparcos, and other catalogues. From this, they calculate 
the open cluster orbits and identified four open cluster com- 
plexes, the members of which share kinetic behaviour and 
age. For example, the youngest complex (age <79 Myr) is 
identified in the Gould Belt. 

Frinchaboy & Majewski (2008) have made a study of 
open cluster 3-d kinematics. The first phase, reported in 
their paper, starts with creating a uniform set of kinematic 
parameters. For this, they measured the radial velocity of 
over 3,000 stars in open clusters and concentrated on those 
with Tycho2 proper motion values. This amounts to 66 open 
clusters, with membership having been determined for ra- 
dial velocity, proper motion, and spatial distribution. Appli- 
cations with orbits are due to follow. 

Finally, and of most relevance to our work, is the study 
of Wu et al. (2009). They use the previous release (2.9) 
of DAML, dated 2008 April 13. Their data analysis ad- 
dresses distances, proper motion, and radial velocity. They 
adopt the DAML distances, but use the proper motion 
data from Kharchenko's ASCC-2.5 catalogue (Kharchenko, 
2001). They find the effect of the re-reduction of Hippar- 
cos data by van Leeuwen (2007) to be less than the error 
values in DAML. Regarding radial velocity, they substitute 
more recent values of Frinchaboy & Majewski (2008), and 
of Mermilliod et al. (2008) into their data set. They use 
three different axisymmetric Galactic potentials and calcu- 
late the orbits of over 400 open clusters for 5Gyr back in 
time. The main conclusions related to orbits are that the 
velocity dispersion of open clusters increases with age, show- 
ing disk heating. They also find that the orbit eccentricity 
and maximum height above the Galactic plane are affected 
most by uncertainties in distance. Finally, they compare or- 
bital eccentricity of globular clusters, disk giants, disk F and 
G dwarfs, and open clusters, and deduce that about 4% of 
open clusters are thick disk objects. 

This work extends the above by concentrating on un- 
usual clusters, rather than the bulk of the population. In this 
paper we study the orbits of the DAML open clusters within 
the Galactic potential. Section 2 presents the available data 
on open clusters, and describes the method of orbit inte- 



gration. In Section 3, we present the results of a catalogue 
data validation, based mainly on orbit morphology, and self- 
consistency tests. In Section 4, we carry out an analysis of 
cluster orbits, and in Section 5 we identify clusters that may 
not be of purely thin disk origin, but instead due to events 
such as the impact on the Galactic disk of a high velocity 
cloud, of a globular cluster, or a merger with a dwarf galaxy. 



2 OPEN CLUSTERS AND ORBIT 
CALCULATIONS 

2.1 Catalogue data 

DAML was developed from an original source by Mermilliod 
(1995), now called WEBDA Q which is still being maintained 
and updated, and which, combined with SIMBAD's proper 
motion data, provides the necessary data for kinematic stud- 
ies of open clusters. These proper motion data are from 
Loktin & Beshenov (2003), and Kharchenko et al. (2005). 
WEBDA provides no errors on distances, but the SIMBAD 
proper motion data have errors attributed to them. 

DAML is a compilation of the data in the literature, of 
which one advantage is a single file with all the key data 
on the open clusters, with the latest updated version 2. 
released on 2009 February 17, which we use throughout this 
work. The DAML catalogue contains entries for 1787 open 
clusters. The data fields in the catalogue include position on 
the sky, distance, proper motion, radial velocity, metallicity, 
size, colour excess, Trumpler type, numbers of stars used in 
determining proper motion, metallicity and radial velocity, 
and where relevant, the error on these quantities. The ex- 
tent of kinematic information available is that, of these 1787 
clusters, 1114 clusters have at least an entry for the distance 
to the Sun, and of these in turn, 481 have the data on dis- 
tance, proper motion, and radial velocity, with 439 of these 
also having the age. Information on metallicity is present for 
111 of the set of 439. 

Of the set of 439, 86% have an age below 0.5Gyr. 
The oldest reported in DAML, Berkeleyl7, is ~ lOGyr, the 
youngest are ~lMyr. The metallicity [Fe/H] range is from 
-0.835 to 0.27 dex, with 20% of values less than -0.2 dex. 

Both WEBDA and DAML, through studies such as that 
of Kharchenko et al. (2005), rely on stellar observational 
proper motion data in surveys, such as Hipparcos (Perryman 
et al. 1997), Tycho2 (H0g et al. 2000), UCAC2 (Zacharias et 
al., 2004), or 4M (Frinchaboy & Majewski 2008). Caution is 
needed with UCAC2 and 4M, because, according to Frinch- 
aboy & Majewski (2008), there may be systematic errors in 
the proper motion values. 

Determination of errors in basic kinematic data is es- 
sential in assessing the strength of any conclusion resulting 
from orbit calculations. Paunzen & Netopil (2006) investi- 
gated the accuracy of other relevant open cluster data such 
as age and distance. They used eight papers on data compi- 
lations, and WEBDA's references, checking that no piece of 
data is used more than once. This then represented data for 
395 clusters, amounting to 6437 individual pieces of data. 
They determined the distributions of the absolute errors in 

1 http://www.univie.ac.at /webda/navigation.html 

2 www.astro.iag.usp.br/~wilton 
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age and distance via a statistical analysis of the sample. 
This shows the peak distribution for age and distance er- 
rors to lie in the respective intervals 40-60%, and 5-10%. 
Paunzen & Netopil (2006) also compared their result with 
the averages given in the 2005 October version of DAML, 
and conclude that if using the parameters from this cata- 
logue, the expected errors are comparable to those obtained 
independently from the literature. This provides support for 
using directly the uncertainty values from the catalogue, a 
view shared by Lepine et al. (2008). 

We recognise that the open cluster membership in 
DAML is in no way volume, or age-limited, and no attempt 
is being made here to study the overall Galactic open clus- 
ter system. We also recognise that the body of data on open 
clusters is changing, with new data being published between 
revisions by Dias and collaborators. Some more recent high- 
quality data, for example on radial velocity measurements 
(e.g. Mermilliod, Mayor & Udry, 2008, and Frinchaboy & 
Majewski, 2008), have not yet been taken into account in 
the latest (2.10) revision of the catalogue. This is an impor- 
tant area of improvement, as Bovy et al.(2009) point out 
that the impact of a single star's radial velocity on that 
of a moving group can be significant. Recognising that in- 
cremental improvements to the data on open clusters are 
taking place, nevertheless, our work aims to reach conclu- 
sions on the set of hitherto observed Galactic open clusters 
for a clearly defined sample as provided in DAML, and we 
therefore do not include these isolated changes. 

2.2 Orbit calculation 

To calculate the (proper, rather than osculating) orbits of 
the clusters, we use a Galactic potential proposed by Fell- 
hauer et al. (2006) for this type of calculation, and assume 
it to be fixed in time. This potential has also been used by 
several authors, for example Law et al. (2005), Fellhauer et 
al. (2007), Vande Putte & Cropper (2009), and Vande Putte, 
Cropper & Ferreras (2009) to track a variety of objects in 
the Galaxy. 

Briefly, the potential is axisymmetric, and therefore uses 
cylindrical coordinates R, </>, z where R is the distance to the 
z axis, 4> is the azimuthal angle between the x axis and the 
projection of the position vector on the xy plane, and z the 
distance of the particle above the xy plane. 

This potential has three parts, representing a bulge, 
halo, and disk: 

$ = $ b + $> h + $ d (1) 

The bulge potential is given by a Hernquist (1990) ex- 
pression: 

GM b 



*f,(r) = - 



r + a 



(2) 



where M b = 3.4 x 10 10 M Q , and a = 0.7 kpc, and where r is 
the Galactocentric radius (r 2 = R 2 + z 2 ). A halo contributes 
to the potential in a logarithmic form: 

2 

= ^ ln(i? 2 + z 2 q- 2 + d 2 ) (3) 

with vo = 186km s _1 , q is a flattening parameter, taken as 
unity here, and d — 12kpc. A Miyamoto-Nagai expression 
(1975) represents the disk contribution: 




Figure 1. Typical crown orbit, for open cluster NGC5316, cal- 
culated for 5Gyr only, for clarity 
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{K 2 + [a + (z 2 + 6 2 ) 1 /2] 2 } 



21 1/2 



(4) 



where M d = 10 11 Mo, a = 6.5kpc, b = 0.26kpc. For this 
potential, the distance from Sun to Galactic centre is set at 
8 kpc, with a local circular velocity of 220km s _1 . 

All calculations are carried out with a fixed Cartesian 
Galactocentric right-handed coordinate system. The x-axis 
points in the direction of motion of the Local Standard of 
Rest (LSR).The y-axis points from the Galactic centre to 
the Sun, so that the z-axis points to the North Galactic 
Pole (NGP). The current cluster positions (x, y, z) used 
in the calculations are obtained from the cluster heliocen- 
tric distance, and its celestial coordinates. The velocities of 
a cluster with respect to the LSR are determined from its 
proper motion and radial velocity, and from the Solar mo- 
tion (U Q = 10±0.4, V Q =5.2±0.6, W =7.2±O.4 km s _1 from 
Dehnen & Binney (1998). Within the LSR's own coordi- 
nate system, the U axis points towards the Galactic cen- 
tre, the V axis points in the direction of motion of the 
LSR, and the W axis points to the NGP. The transfor- 
mation from radial velocities and proper motions to mo- 
tion in the LSR (U, V, W), uses the formalism of Johnson 
& Soderblom (1987). The elements of the transformation 
matrix of Johnson & Soderblom have been updated to the 
International Celestial Reference System (ICRS), based on 
RA(NGP)= 192.85948°, Dec(NGP)=+27.12825°, and the 
ascending node of the Galactic plane on the equator, Iq= 
32.93192° (ESA, 1997). The transformation from velocities 
in the LSR to velocities in our fixed Galactocentric sys- 
tem (dx/dt= V+ 220 km s _1 , dy/dt = -U, dz/dt = W) 
uses an LSR velocity of 220±15 km s~\ The origin of the 
LSR is placed at 8kpc from the Galactic centre (Binney & 
Tremaine, 2008). 

The relevant equations of motion are integrated using 
a fourth-order Runge-Kutta procedure implemented in a 
FORTRAN code, and the code tolerance parameters set so 
as to conserve total energy and angular momentum to better 
than 10" 10 . 

We show in Table [T]some examples of the observed data 
relevant to the orbit calculations, whereas Table [5] displays 
the corresponding input data in the case of our Galactocen- 
tric representation. In Table [2] the distance error is set at 
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10%, and the error on the velocities also accounts for the 
error in the Sun's motion in the LSR, and in the circular 
velocity at the Solar radius. All observed data in Table [1] 
are from DAML (version 2.10). 



3 DATA VALIDATION 

In this section we present the results of an initial screening 
of the data relating to those clusters for which sufficient 
information is available for orbit calculations. 

We considered the 481 clusters that have at least dis- 
tance, proper motion, and radial velocity data, and com- 
puted the orbits backwards for up to 15Gyr, visually exam- 
ining each orbit. The value of 15Gyr is longer than the age 
of the universe, and we recognise that the potential will have 
varied during this time, but this arbitrarily long time is used 
to place all cluster orbits on an equal basis. Most clusters 
show quasi-periodic crown orbits, an example of which ap- 
pears in Figure [1] in this case, for NGC5316. There are four 
exceptions: Berkeley 20, Berkeley 29, Berkeley31, and Berke- 
ley33 do not produce typical crown orbits. Figure [5] displays 
their orbits after 15Gyr back in time and Figure [3] shows sky 
images and positions of these four clusters. They also share 
the distinction of being the clusters reaching the highest ver- 
tical distance from the disk plane (called z max ), calculated 
over the orbit described in 15 Gyr (Table |3J. Furthermore, 
these clusters' absolute value of nominal z at birth is high- 
est of all, with z at birth being -3.8 kpc for Berkeley20, 13.9 
kpc for Berkeley29, 12.9 kpc for Berkeley31, and 6.7kpc for 
Berkeley33. The four clusters have remarkable properties: 
they are among the 10% oldest for those 439 with known 
age, are among the seven most distant of our sample of 481, 
and have among the twenty lowest metallicity values in our 
sample of 111. Clearly, the four clusters exhibit an interest- 
ing kinematic behaviour, meriting closer examination. We 
summarise their properties in Table [3] 

3.1 Sense of Rotation 

We note that Berkeley 29, one of the four clusters singled 
out in the pre-screening, counter-rotates with respect to the 
Galactic rotation: it is the only one of the 481 in the sample 
that does. Uncertainties in the values of the Sun's motion in 
the LSR, velocity of the LSR, distances, radial velocities, and 
proper motions affect the cluster orbit (including its sense 
of rotation). To examine the effect of these uncertainties on 
the sense of rotation of Berkeley29, we assumed that the un- 
certainties for Solar motion, LSR motion, distances, radial 
velocities and proper motions are all normally distributed, 

3 Since this work was submitted, new values for the Solar mo- 
tion have been proposed (Schonrich, Binney & Dehnen, 2010): 
(U,V,W) e = (ll.lj;g;™,12.24±S;g,7.35±S;|;)kni s~\ with ad- 
ditional systematic uncertainties ~ (1,2,0.5) km s" 1 . We note 
that, Vq is 7 km s" 1 larger than previously estimated (see above), 
and this will have the greatest effect on dx/dt in Table [2] How- 
ever, given the magnitude of dx/dt and the size of the uncertainty 
in the velocity of the LSR (~15 km s _1 ), little effect is found on 
the nature of the orbits, or on the orbital parameters z max and 
r\ used in Table [4] (differences for the top seven candidates for 
external origin of a few percent of the uncertainties). 



with standard deviations equal to the quoted experimental 
errors from the Dias database. We then performed a Monte 
Carlo simulation for Berkeley29, with 1,000 trials. The result 
is that for 84% of cases, the rotation is counter to that of 
the Galactic disk. We note that Frinchaboy (2009) found a 
similar orbit to that of Figure [21 also counter-rotating, using 
the kinematic information from the DAML catalogue. While 
stars from a cluster originating in a cannibalised small com- 
panion could have counter-rotating orbit, counter rotation 
of stars is generally unusual, being observed in less than 1% 
of disk galaxies (Sparke & Gallagher, 2000). 

3.2 Coordinate Differences 

We notice that there is a small difference in coordinates 
between SIMBAD and the DAML catalogue for Berke- 
ley20: the SIMBAD coordinates (RA; Dec, J2000) are: a = 
05 h 33 m 00 s ; S = +00°13'00", whereas in DAML these val- 
ues are a = 05 h 32 m 37 s ; S = +00°11'18", which represents 
a separation of 5.99 arc min. This, however, proves to be 
insufficient to produce the strange orbit in Figure [2] 

We notice that there is also a small difference in co- 
ordinates between SIMBAD and the DAML catalogue for 
Berkeley29: the SIMBAD coordinates (RA; Dec, J2000) are: 
a = 06 h 53 m 04:2; S = +16°55'39", whereas in DAML these 
values are a = 06 h 53 m 18 s ; 5 = +16°55'00", which repre- 
sents a separation of 3.35 arc min. 

For Berkeley31 and 33, the SIMBAD and DAML coor- 
dinates are identical. Hence the cause of the odd orbits does 
not lie with the small coordinate differences. This leaves the 
distance and proper motion as possible sources of error. 

3.3 Proper Motions 

The proper motion data in DAML are based mainly on the 
work of Dias et al. (2006), Kharchenko et al. (2003), and 
Loktin (2003), using the UCAC2 catalogue (Zacharias et al., 
2004), the Tycho2 catalogue (H0g et al, 2000), or the Hippar- 
cos catalogue (Perryman et al, 1997), with cluster positions 
being identified by the centroid of a set of coherently-moving 
stars. In the case of Berkeley20, 29, 31, 33, the proper mo- 
tions are determined (respectively) by Dias et al. (2006) from 
UCAC2, by Loktin & Beshenov (2003) from Tycho2 data, 
by Dias et al (2002a) from Tycho2, and by Dias et al. (2006) 
from UCAC2. 

Close scrutiny of the data relating to the individual 
stars used in DAML to determine the proper motion of 
Berkeley20 reveals that their position is ~ 74° away from 
the cluster. This mismatch, probably due to a clerical error 
in the catalogue, explains the strange orbit. 

3.4 Distance 

There is some weak evidence for the distance to these par- 
ticular clusters of "proper motion" stars being smaller than 
given in the catalogue for Berkeley29, 31, and 33. Firstly, 
some estimate of the distance from Berkeley29 to the set of 
stars used for proper motion determination is possible using 
main sequence fitting. FigureUshows the B & V magnitudes 
for these stars, based on the GSC2 catalogue (Lasker et al, 
2008). Most stars are on the main sequence, and based on 
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Table 1. Observed data relevant to orbit calculations for the first ten open clusters in DAML with sufficient data to calculate orbits, 
with data extracted from DAML. A full version of the table for all 481 clusters for which an orbit calculation is possible, appears in 
the electronic version of the paper. 
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Berkeley 59 


00 : 


02 


: 14 


+67 ; 


: 25 : 


00 


1.000 ± 0.200 


-12.5 ± 7.1 


-2.11 ± 0.81 


-1.20 + 0.75 


Blanco 1 


00 : 


04 


: 07 


-29 : 


50 : 


00 


0.269 + 0.054 


4.1 ± 1.4 


20.17 + 0.51 


3.00 + 0.51 


Alessi 20 


00 : 


09 


: 23 


+58 : 


39 : 


57 


0.450 ± 0.090 


-11.5 ± 0.01 


8.73 + 0.53 


-3.11 + 0.53 


ASCC 1 


00 : 


09 


: 36 


+62 : 


40 : 


48 


4.000 ± 0.800 


-69.7 + 4.7 


-2.07 + 0.72 


0.46 ± 0.57 


Mayer 1 


00 : 


21 


: 54 


+61 : 


45 : 


00 


1.429 + 0.286 


-20.9 ± 2.0 


-4.46 ± 1.13 


-6.66 + 0.94 


NGC 129 


00 : 


30 


: 00 


+60 : 


13 : 


06 


1.625 + 0.325 


-37.4 + 0.5 


-1.06 ± 0.94 


1.60 + 0.94 


ASCC 3 


00 : 


31 


: 09 


+55 : 


16 : 


48 


1.700 + 0.340 


-37.0 + 0.0 


-1.92 ± 0.61 


-1.25 + 0.59 


NGC 225 


00 : 


43 


: 39 


+61 : 


46 : 


30 


0.657 + 0.131 


-28.0 ± 0.0 


-4.95 ± 0.76 


-0.50 + 0.76 


NGC 188 


00 : 


47 


: 28 


+85 : 
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18 


2.047 + 0.409 


-45.0 + 10.0 


-1.48 ± 1.25 


-0.56 ± 1.24 


IC 1590 


00 : 


52 


: 49 


+56 : 


37 : 


42 


2.940 ± 0.588 


-32.5 ± 6.4 


-1.36 ± 0.23 


-1.34 + 0.83 



Table 2. Present positions and velocities of open clusters in Table[T] as calculated for use in the orbit code. A full version of the table 
for all 481 clusters for which an orbit calculation is possible, appears in the electronic version of the paper. 
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X 




y 

(kpc) 


z 


dx/dt 


dy/dt 
(km s -1 ) 


dz/dt 


Berkeley59 


0.878 


± 


0.176 


8.471 


± 


0.094 


0.087 + 0.017 


219.7 


± 


16 


-25.3 + 5.1 


2.4 + 3.7 


Blancol 


0.013 


± 


0.003 


7.952 


± 


0.010 


-0.264 ± 0.053 


217.2 


± 


15 


13.4 + 4.9 


-1.7+1.7 


Alessi20 


0.398 


± 


0.080 


8.207 


± 


0.041 


-0.030 ± 0.006 


206.5 


± 


15 


-0.2 + 3.2 


-1.6 + 2.3 


ASCC1 


3.527 


± 


0.705 


9.887 


± 


0.377 


0.014 + 0.003 


181.4 


± 


17 


-75.8 + 13.9 


21.9 ± 11.3 


Mayerl 


1.244 


± 


0.249 


8.702 


± 


0.140 


-0.023 ± 0.005 


223.8 


± 


16 


-51.2 + 9.1 


-33.8 + 10.5 


NGC129 


1.402 


± 


0.280 


8.818 


± 


0.164 


-0.072 + 0.014 


197.1 


± 


15 


-34.7 + 6.4 


21.8 + 7.7 


ASCC3 


1.459 


± 


0.292 


8.843 


± 


0.169 


-0.221 + 0.044 


200.6 


± 


15 


-43.0 + 5.2 


3.2 + 5.0 


NGC225 


0.557 


± 


0.111 


8.348 


± 


0.070 


-0.012 + 0.002 


209.7 


± 


15 


-38.0 + 3.3 


6.6 + 2.4 


NGC188 


1.590 


± 


0.318 


9.027 


± 


0.205 


0.780 + 0.156 


199.8 


± 


19 


-43.6 + 11.8 


-14.8 ± 11.8 


IC1590 


2.448 


± 


0.490 


9.597 


± 


0.319 


-0.320 ± 0.064 


206.8 


± 


16 


-44.6 ± 5.6 


-8.0+12.1 



the position of the stars in the CM diagram, we believe that 
the stars along the oblique could be cluster stars, with the 
rest being field stars. This oblique is almost parallel to the 
Hyades sequence, so this allows us to use the main sequence 
fitting method to find the distance to the cluster, based on 
the known distance to the Hyades cluster (see for example 
Karttunen et al., 1996). This uses the relation 

V = Mv + 51og(^-) +A V (5) 

where V refers to the cluster at a given B — V, My to the 
Hyades cluster at the same B — V, with Av determined from 
the NED databas^l which incorporates the data of Schlegel 
et al. (1998). We estimate the distance to the cluster con- 
taining the stars used in proper motion determination to be 
~lkpc rather than the DAML value of 15kpc (Table [3]). 

Frinchaboy (2009) further concluded that the UCAC 
stars used to determine the proper motion of Berkeley29 are 
too bright to be part of that distant system. A similar ar- 
gument to that made for Berkeley 29 produces a distance of 
~lkpc for the stars used in proper motion determination of 

4 http://nedwww.ipac.caltech.edu/ 



Berkeley31. No equivalent value is available for Berkeley 33, 
as no V magnitudes are given in GSC2 (Lasker et al, 2008) 
for the stars selected in the proper motion determination. 

We then checked that reducing the distance to Berke- 
ley29 to 1 kpc produces a crown orbit for the DAML proper 
motion and radial velocity, and we find that the cluster now 
rotates in the Galactic sense. Likewise, reducing the distance 
to Berkeley31, and 33 to lkpc also results in crown orbits. 
This suggests that the proper motion values in DAML might 
refer to an association close by, distinct from Berkeley29 (or 
Berkeley31, or Berkeley33). 

A second piece of evidence for a mismatch of distance 
and proper motion is provided by an estimate of proper mo- 
tion for an object at the same distance and direction as the 
cluster, but in the disk, and assuming a constant circular ve- 
locity. This is established using a formalism in Sparke and 
Gallagher (2000) and yields values of 0.01, 0.05, and 0.25 
mas y _1 for Berkeley29, 31, 33, which is 1-2 orders of mag- 
nitude less than the values given in the catalogue (see Table 

As a third piece of evidence, we note that Berkeley20, 
29, 31, 33 are the most, second most, third most, and seventh 
most distant clusters from the Sun (Table [3]). At these large 
distances, the cluster members may be too faint for the pro- 
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Figure 2. Orbits for 15 Gyr back in time 



Table 3. The four clusters that do not describe crown orbits. Columns 2 to 6 arc from DAML 



Cluster 


fi a cos5 


US 


d & 


Age 


[Fe/H] Solar 


z max 




(mas yr _1 ) 


(mas yr _1 ) 


(kpc) 


(Gyr) 


(dcx) 


(kpc) 


Berkeley20 


1.51±0.81 


-4.11±0.82 


8.4±1.68 


6.02±1.51 


-0.61±0.14 


7.63±4.11 


Berkeley29 


-0.14±0.8 


-4.75±0.58 


14.87±2.97 


1.06±0.26 


-0.31±0.03 


19.9±18.68 


Berkeley31 


-4.3±0.52 


-3.97±0.52 


8.27±1.65 


2.05±0.51 


-0.40±0.16 


15.2±10.76 


Berkeley 33 


-5.73±0.96 


3.71±0.96 


6±1.2 


0.79±0.19 


-0.26±0.05 


7.97±11.21 



Note that Carraro et al. (2007) give a very different value for the age of Berkeley29, namely 4.5Gyr. 



cess of proper motion star selection to have detected them, 
and here we note that the limiting magnitude of UCAC2 
(or even UCAC-3) is 16 (Zacharias et al., 2004, Zacharias et 
al., 2010, and Finch, Zacharias & Wycoff, 2010). This would 
then favour stars closer by. 

In principle, proper motion information for fainter stars 
are available in the USNO-B1.0 catalogue (Monet et al. 
2003), and the recently released PPMXL catalogue (Roser, 
Demleitner, Schilbach, 2010). However, their proper motion 
accuracy is not better than ~ 4 mas yr _1 , which is larger 
than the expected proper motions for the cluster stars (see 
above) . 



3.5 Evidence from Non-kinematic Indicators 

The unusually old ages and low metallicities of Berkeley 20, 
29, 31, and 33 have been highlighted in the literature, for 
example in Yong et al., (2005), Carraro & Bensby (2009), 
MacMinn et al. (1994), and Carraro et al. (2007). The values 
are included in Table [3] 



3.6 Adopted Sample 

Our orbit screening has identified possible anomalies in the 
data for Berkeley 20, 29, 31 and 33. Further analysis of the 
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Figure 3. Positions of the four clusters Berkeley20, 29, 31, 33 according to SIMBAD (red dot), and DAML (yellow dot). SIMBAD and 
DAML positions coincide for both Berkcley31 and Berkcley33, so only a red dot is shown. The red open circles have a radius of 25" and 
are centered on the 10 stars that have the highest probabilty of cluster membership. There are no open circles for Berkeley20, as the 
proper motion stars lie 74° off the nominal position of the cluster. The figures are ~19' wide and ~13' high, and are sourced from the 
DSS at http:archive.stsci.edu/cgi-bin/dss_form 



input data for the orbit analysis has shown that the data 
may be in error. We classify Berkeley20, 29, 31, 33 as out- 
liers on-hold because of their orbit appearance, their high 
proper motion in relation to their distance in the Galac- 
tic potential, and because their large distance will bias the 
process of selecting proper motion stars against the fainter 
clusters stars, in favour of closer associations of field stars. 
For Berkeley29 and 31, proper motion stars could be fore- 
ground objects much closer to us. We plan to revisit these 
four clusters when newer data on proper motion have been 
carefully analysed and cluster proper motion measured. 

Having found the small mismatch in coordinates be- 
tween SIMBAD and DAML for Berkeley29, we then carried 
out a comparison between Simbad and DAML coordinates 
for the 368 out of 481 clusters in DAML clusters that are 
resolved by SIMBAD. The comparison shows 3 cases where 
the difference between Simbad and DAML is in excess of 1°: 
Markarian38 (113°), Markarian6 (27°), and ASCC14 (54°). 
We remove these three objects from further consideration, 
particularly as the SIMBAD database lists the Markarian 
objects as galaxies, and the third as a single, high-proper 
motion star. 

Based on the observed large difference between the clus- 
ter position in DAML, and the centroid of the set of stars 
used to determine proper motion for Berkeley20(~ 74°), we 
undertook a comparison of the celestial coordinates of the 



centroid of stars used to determine proper motion, versus the 
coordinates in the main DAML database. To the extent that 
the positions of the stars used in proper motion studies are 
given in the Tycho2 (H0g et al, 2000), TJCAC2 (Zacharias et 
al, 2004), and Hipparcos (Perryman et al, 1997) catalogues, 
we were able to effect this comparison for 252 out of the 481 
clusters. We found that apart from Berkeley20 already noted 
above, there is only one other cluster where the difference 
exceeds 1°: Bochum4 (130°). This cluster is also removed 
from further consideration. 

Finally, we remove from consideration clusters for which 
less than five stars are used to determine proper motion. 
This involves rejecting 63 clusters from our sample of objects 
with sufficient data to calculate orbits, leaving 410 clusters. 
The resulting subset of clusters with also an age contains 
369 members, of which 96 also have metallicity data. 



4 ORBIT ANALYSIS 

We now look in more detail at orbits and metallicity in our 
cluster population. 

Most clusters have orbits close to the Galactic plane. 
We expect such behaviour from clusters born from disk gas 
condensation, with their velocity lying in the plane of the 
disk, and aligned with the local circular velocity. 
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Figure 4. CM diagram for stars used in proper motion determi- 
nation of Berkeley29 in DAML. The oblique runs through pre- 
sumed cluster stars (see section 3.4) 



A parameter of interest in this respect is z m ax , the max- 
imum distance of the orbit from the Galactic plane. This 
shows whether the orbits could be associated with the thin 
or thick disk, or beyond. We can thus view z max as a diag- 
nostic of non-alignment with the Galactic plane velocities. 

We can also use the radial quantity 

Rmax Rmin i r .\ 

71 = Tr T~r — W9 (6) 

where the variable R is the distance between the cluster 
and the Galactic centre, projected onto the disk. The param- 
eter rj (twice the eccentricity for an elliptical orbit) will be 
small if the cluster velocity is aligned with the local circular 
velocity. It is therefore also a diagnostic for non-alignment 
with the circular velocity. 

Figure [S] shows the cumulative distribution of rj and 
z m „ values, illustrating the fact that most open clusters 
have low z ma x and rj, as would be expected for an origin 
in the disk. There are clear inflexions in the distribution, 
at z max =0.35 kpc (90% of clusters) and ?7=0.24 (80% of 
clusters), indicating a possible transition to a different pop- 
ulation. In order to select a sample of the most extreme 
examples well beyond these inflexion points, we set the se- 
lection thresholds for segregating our results at z max =0-9 
and 77=0.5. These values correspond to z„ lax —3 times the 
scale height of the thin disk (Binney & Merrifield, 1998), 
and rj when the orbit thickness in the disk plane is half 
the average distance between orbit and Galactic centre. We 
consider that values above these thresholds will indicate an 
unusual origin. 

The metallicity is clearly also a diagnostic of past his- 
tory, and we also define a metallicity of -0.2 above which 
we find most clusters, to be a threshold below which we 
presume an unusual origin for a disk population. 

We plot the various combinations of z max , rj, and metal- 
licity in Figures [6] and [7| Note that the metallicity error bars 
in Figures [5] and [7] represent the errors given in DAML, and 
in the case of z max and rj are the standard deviation obtained 
from a Monte Carlo calculation carried out for 1000 realisa- 
tions per cluster for the 410 clusters that satisfy the initial 
screening, using the error estimates in DAML as standard 
deviations of normal distributions, and a past orbit integra- 



tion time of 15Gyr. There is a correlation between all three 
parameters, z max , rj, and metallicity. 

Table 2] shows the main data for clusters that meet at 
least one of the criteria for unusual origirjf]. We see that of 
the 410 clusters initially considered, 35 satisfy at least one 
criterion. Four satisfy all three criteria for unusual origin: 
Berkeley21, 32, 99, and Melotte66, and here the z max versus 
rj linear correlation coefficient is 0.97. These clusters are aged 
2.19, 3.38, 3.16, and 2.78 Gyr, respectively, which places 
them amongst the ten oldest clusters among the set of 369 
for which ages are available in DAML. Berkeley21 has the 
highest Zma X and rj, and its orbit for 15 Gyr is thick and 
wide. Thirteen others have data on metallicity, and of these, 
four meet two of the criteria, whereas nine meet only one 
criterion for unusual origin. We looked for information on 
cluster masses, but found only eleven clusters for which the 
mass is available (Durgapal & Pandey, 2001), of which only 
two values are for our set of clusters in Table 2] At least for 
these, the masses are not atypical in any degree. In DAML, 
we find information on actual diameters in pc, but for the 
clusters in Table U we find no correlation between these 
diameters with age or metallicity. 

Finally, as support evidence, we have calculated the po- 
sition of the birthplace of clusters which meet all three cri- 
teria using a Monte Carlo method with 2000 realisations. 
Figure [8] shows the results in terms of z and 7? at birth for 
the four prime candidates for unusual origin. For Berkeley21, 
there is a 92% chance of birth occurring at z ^ 0.3 kpc (one 
scale height). For Berkeley32, this figure is 72%. For Berke- 
ley99, there is a 90% chance of birth beyond 0.3 kpc from 
the disk plane, increasing to 94% if we also require that R be 
less than the visible disk diameter of 15kpc. For Melotte66, 
these figures are 88% and 92% respectively. The nominal val- 
ues of z at birth are for Berkeley21 2.5kpc, for Berkeley32 
0.7kpc, for Berkeley99 -1.8kpc, and for Melotte66 -1.3kpc; 
these are among the ten highest distances from the plane 
calculated for those in our sample of clusters that also have 
age data. In addition we carried out a similar calculation for 
NGC5316 (Figure [T]), and find that all values of z at birth 
lie within the thin disk with a FWHM of 0.4kpc. 



5 CLUSTER ORIGIN 

We now examine the nature of the unusual origin behind 
the clusters in Table [4] by considering the mechanisms that 
could lead to the formation of open clusters. We see that 
most clusters in Tabled and certainly those where at least 
two criteria are met, exhibit relatively high values of z max - 
For high-altitude clusters such as these, two formation mech- 
anisms have been suggested. Martos, Allen, Franco & Kurtz 
(1999) modelled the gas response to the spiral arm density 
wave and find that the shock sends gas from the inter-arm 
region to high altitude, followed by star formation, albeit 
with a low efficiency. The other mechanism, proposed by de 

5 As stated above, Wu et al., (2009) use an earlier version of the 
Dias catalogue complemented by some third-party proper motion 
data; they also use a different expression for the gravitational 
potential to the one used here (Allen & Santillan, 1991) so some 
differences in z ma x and rj can be expected, but we find these to 
be well within the uncertainties on our values. 
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Figure 5. Cumulative probability distribution of z ma x, left panel, with threshold value of 0.9 shown by dashed line, and cumulative 
probability distribution of r\ right panel, with threshold value of 0.5 represented by dashed line 
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n 



Figure 6. Relation between z m ax and r\. The green lines at z maa; =0.9 and r)=0.5 represent the thresholds beyond which unusual origin 
is a possibility. The darker grey area represents the zone of foremost interest for clusters of unusual origin, with two criteria for unusual 
origin satisfied, the lighter grey corresponds to areas for which only one of the criteria is satisfied 



la Fuente Marcos & de la Fuente Marcos (2008) has pre- 
existing clusters, preferably globular, interacting with high 
altitude gas clouds, creating enhanced turbulence leading to 
star formation in sufficiently large clouds. They further es- 
timate that about 3,000 clusters could have formed in this 
way over the life of the Galaxy. 

The above mechanisms could explain most clusters in 
Table [4] However, we offer further alternative mechanisms 



for consideration. Firstly, as discussed by Friel (1995), im- 
pact of high velocity clouds on the disk can lead to star for- 
mation, a mechanism studied by Comeron & Torra (1992), 
with the cluster retaining some kinetic memory of the event, 
for example high z (Danly, 1992) and possibly eccentric or- 
bits. In addition, Wakker et al. (1999) find that such High 
Velocity Clouds have low metallicities (~0.1 Solar), with 
some of them falling into the Galaxy. These point to the 
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Figure 7. Relation between z m ax, and metallicity (top), and between r) and metallicity (bottom) The green lines represent the Zmax 
and r) thresholds beyond which unusual origin is a possibility. The green line at metallicity —0.2 is the threshold below which unusual 
origin is a possibility. The darker grey area represents the zone of foremost interest for clusters of unusual origin, with two criteria for 
unusual origin satisfied, the lighter grey corresponds to areas for which only one of the criteria is satisfied 



origins of clusters with high z ma x and rj and low metallic- 
ity being formed as a result of an agent of extra-Galactic 
origin: in our case, these would be the top four clusters in 
TableU and might also be the case for NGC2158, 2420, 7789, 
and IC1311. Secondly, globular clusters impacting the disk 
can cause disk gas compression, either due to gravitational 
focussing (initially discussed in Wallin, Higdon & Staveley- 
Smith 1996, and later in Vande Putte & Cropper, 2009), 
or shock wave formation (Levy, 2000). This would produce 
clusters with a metallicity representative of the disk, and 
the parameters z max and 77 would be the indicators of un- 
usual origin; in Table|4]these might be NGC6791, 1817, and 
7044. Next, the accretion of satellite galaxies could result 
in donations of ready-made open clusters with a metallic- 



ity higher than that of the high velocity clouds. Here too, 
z mol and 77 would be diagnostics, because of the cluster re- 
taining kinematic memory of the donor galaxy. Metallicity 
could also be a diagnostic, but is likely to be higher than 
that of the high velocity clouds. Finally, clusters could be 
formed as the result of a merging galaxy, such as within the 
dwarf Sagittarius or the Magellanic steams. Here we would 
expect again high z ma x and 77 and perhaps a higher or sim- 
ilar metallicity. Some unspecified form of merger could be 
the agent responsible for the last nine clusters in Table 2] 
Clusters originating from these last two mechanisms should 
be related at least in terms of their integrals of motion to 
Galactic streams and/or to known HI features such as the 
Magellanic stream. 
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From the literature we find the following non-kinematic 
studies on the origins of these four clusters and for those 
we have placed "on hold". Concerning our four prime can- 
didates for unusual origin identified by orbit analysis, we 
find that Yong et al. (2005) are unable to reach a conclu- 
sion for Berkeley21. There are no further relevant studies for 
Berkeley32, 99, and Melotte66. In the case of the clusters 
"on hold", Yong et al. (2005) using high resolution spec- 
troscopy, provide arguments based on age and metallicity 
to conclude that Berkeley29 could be due to a merger, but 
not Berkeley31. They are unable to reach a conclusion for 
Berkeley20 or 21. Carraro & Bensby (2009) reach the same 
conclusion in respect of Berkeley 29, noting that it lies in 
the dSgr stream. Carraro et al. (2007) conclude that Berke- 
ley33 is not of extra-Galactic origin, but instead is a genuine 
Galactic cluster. 

While some small fraction of currently high Galactic 
altitude clusters could have been born at low altitudes, they 
would be expected to be of at least moderate metallicity. For 
our four prime candidates this is not the case, supporting 
their origin high above the Galactic plane. 



6 CONCLUSION 

The DAML catalogue of 1787 open clusters provides a large 
database to study these clusters. In particular, it contains 
sufficient data to calculate past orbits in as many as 481 
cases. This allowed us to undertake a systematic analysis of 
these open clusters by examining their orbit morphology. 

We found, with few exceptions, that the orbit type, 
when run for a sufficiently long time, is a crown orbit. This 
type of orbit is readily characterised by two parameters: 
Zmax, the largest separation between cluster and disk in ab- 
solute value, and r], a measure of orbit eccentricity. This has 
first of all revealed four clusters whose proper motion data 
in the catalogue are open to question, as they may relate in- 
stead to foreground associations: they don't exhibit a crown 
orbit, and one of them counter rotates. These are Berke- 
ley 20, 29, 31, and 33, four otherwise interesting clusters with 
low metallicity and large distance from the Sun, but for 
which we found a mismatch between positional and kine- 
matic information in DAML. We temporarily place these 
"on hold" for further study when more data on proper mo- 
tion become available, as it is conceivable that our analysis 
has identified some confusion in the proper motion analysis 
reported in the catalogue. This means that there may be 
four as-yet unidentified associations, generally closer to us 
than the four clusters listed above. 

We then eliminate from the cluster sample those with 
less than five stars (63 clusters) and four other clusters where 
data may be affected by clerical error. We then examine the 
relation between z max , rj, and metallicity in the resulting 
sample of 410. 

We suggest that the higher z max and rj, and the lower 
the metallicity, the chances increase of the cluster being of 
unusual origin. We fix thresholds of >0.9kpc, and>0.5 on 
Zmax, V- If we a l so set a further threshold of <-0.2 in the 
metallicity [Fe/H] dex, we find four clusters breaching all 
three thresholds: Berkeley21, 32, 99, and Melotte66 which 
are the four strongest candidates for unusual origin: we ar- 
gue this could be extra-Galactic (resulting from the impact 



of an external high velocity cloud). An additional 24 clus- 
ters breach at least one of the thresholds, and we identify 
four other possible contenders for generation by HVC impact 
(NGC2158, 2420, 7789, and IC1311). Three clusters may be 
due to globular cluster impact on the disk i.e. internal origin 
(NGC6791, 1817, and 7044). Another four could be due to 
some form of merger. These mechanisms are offered as alter- 
natives to two previously suggested mechanisms for cluster 
formation at high altitude. 

The error values given in the catalogue illustrate an- 
other important point. Such uncertainties are bound to af- 
fect any study of the birth of open clusters, and is an exam- 
ple of how the ESA Gaia mission can help to substantially 
improve our understanding of our Galaxy's formation. 
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Figure 8. Monte Carlo realisations (2,000) for the birthplace of the four prime candidates for origin outside the disk. The horizontal 
lines delimit the thin disk (±0.3 kpc), whereas the vertical line delimits the visible disk (15 kpc). 
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Table 4. Main orbital parameters for clusters that meet at least one of the criteria for possible unusual origin. A full version of the table 
(first 5 columns) for all 481 clusters for which an orbit calculation is possible, irrespective of how many criteria for unusual behaviour 
are met, appears in the electronic version of the paper. 



Name 




V 


[Fe/H] Solar 


Age 


Stars 


Crit. 


Mech. 




(kpc) 




(dex) 


(Gyr) 








Berkeley21 


4.22±2.67 


1.06+0.26 


-0.835+0.1 


2.188+0.547 


32 


z,e,m 


Extra-Gal. 


Berkeley99 


2.20+1.77 


0.78+0.25 


-0.58+0.10 


3.162+0.791 


5 


z,e,m 


Extra-Gal. 


Melotte66 


1.72+1.24 


0.62+0.27 


-0.33+0.03 


2.786+0.697 


10 


z,e,m 


Extra-Gal. 


Berkeley32 


0.99±0.68 


0.61+0.15 


-0.29+0.04 


3.388+0.847 


15 


z,e,m 


Extra-Gal. 


NGC2354 


2.14±1.22 


0.65+0.14 


- 


0.134+0.033 


20 


z,e 




NGC1893 


1.51±1.51 


0.93+0.27 


- 


0.003+0.001 


15 


z,e 




NGC6791 


1.37±1.12 


0.58+0.22 


0.11+0.10 


4.395+1.099 


15 


z,e 


GC disk 




n i fi-t-n 08 


KJ. O J- znu. 6 -L 


9^+0 OQ 


1 0^4+0 9fi4 

J_ . KJO*±^KJ . 


20 


e,m 


J_j A L I a Vjdil . 


in vjvj ( i oy 


90-1-0 0^ 


kj . uo^u. ru 


94+0 OQ 
u. £*±^zkj. uy 


1 41 3+0 3^3 


10 


e,m 


Ija L 1 ct VJJ ell . 


IC1311 


7Q+0 ^0 


KJ . U J- ro 


3+0 1 R 


1 ^8^+0 3QR 


25 




JjA L 1 ct VJJ Cll . 


NOC2420 

IN VJT ^ rr & KJ 


86+0 1 8 


n ^8+0 1 3 

D . tJO_l_L/. -L<J 


_0 38+0 07 

kj. - J o i w . kj i 


i nqc;+o 499 


76 


e,m 


1 iA L 1 cl v_J dl ■ 


IN \Ji V_.' zooo 




n 49+0 1 Q 




U. 1ZUIU.UOU 


8 






RorlcTH 4 

1_) \_ 1 ±\A_ 1 K^ y irt 




n 4^+0 1 4 




1 c;or; i n one 

-L . iJOiJ_l_L/. 3U 


29 






NOC2324 


qi +o 80 


n 47+n 20 

U. rl 1 _I_U . AKJ 


_0 1 7+0 05 


447+0 1 1 2 

U.'±'± 1 _I_U . llii 


10 






r)olid7e25 


7Q-I-1 09 

KJ. 1 L . KJtJ 


i 02+0 24 




006+0 002 


25 






ASCC63 


1 9+0 1 8 


qs+0 09 




01 7+0 004 


10 






Pitjm i c; 1 7 

J. lollllo _i_ 1 


1 2+0 1 

KJ. J_ii_l_W. J-W 


fi7+0 1 2 




01 1 +0 003 


20 






1 llt'H S 




n fi^-r-O 1 9 






12 






ASCC43 


n 1 1 _i_n Q4 

KJ. J- J — (Jrl 


n fiA-Un 1 3 

D . UTt_l_LJ. J-iJ 




i qi +o 048 

VJ. J-is J — 1_ KJ . KJ'-tKJ 


19 






r 1 Vi i in nl^rl h 
±1 UliljJlCI _LU 


n flK-r-D 9^ 


U. U J. IT_U. TO 




U.UUJUJ.UUJ. 


8 






NGC436 


0.89+0.40 


0.59+0.13 


- 


0.084+0.021 


17 


e 




Bochum2 


0.03+0.09 


0.58+0.17 




0.005+0.001 


20 


c 




NGC7044 


0.69+0.84 


0.57+0.19 


-0.16+0.09 


1.901+0.475 


10 


c 


GC disk 


NGC1817 


0.78+0.29 


0.57+0.16 


-0.14+0.05 


0.409+0.102 


10 


e 


GC disk 


Ruprecht47 


0.09+0.15 


0.54+0.13 




0.078+0.019 


7 


e 




ASCC17 


0.24+0.11 


0.53+0.16 




0.013+0.003 


11 


c 




NGC2506 


0.65+0.18 


0.37+0.13 


-0.2+0.02 


1.109+0.277 


20 


in 


Merger 


Kingll 


0.36+0.14 


0.19+0.11 


-0.27+0.15 


1.117+0.279 


10 


in 


Merger 


NGC869 


0.18+0.06 


0.21+0.11 


-0.3 


0.012+0.003 


4 


in 


Merger 


NGC884 


0.20+0.08 


0.11+0.08 


-0.3 


0.013+0.003 


18 


in 


Merger 


NGC6716 


0.14+0.03 


0.09+0.07 


-0.31 


0.091+0.023 


11 


in 


Merger 


NGC2204 


0.78+0.15 


0.41+0.12 


-0.32+0.1 


0.787+0.197 


19 


in 


Merger 


IC2581 


0.28+0.19 


0.12+0.08 


-0.34+0.06 


0.014+0.003 


10 


in 


Merger 


NGC6451 


0.43+0.18 


0.16+0.11 


-0.34 


0.136+0.034 


9 


in 


Merger 


NGC2451B 


0.05+0.02 


0.05+0.08 


-0.45 


0.044+0.011 


5 


in 


Merger 



§ "Stars" represents the number of cluster stars used to determine the proper motion. "Crit" is the criterion for unusual origin, with e 
representing r}, z referring to z max , m being metallicity. A dash indicates that the information is not available. "Mech" refers to the 
formation mechanism, with Extra-Gal meaning an extra-Galactic agent has caused the cluster to form in the Galaxy, GC disk referring 
to a globular cluster impacting gas in the disk, and Merger referring to donation from a infalling galaxy. The small horizontal gap in 
the table separates the clusters for which both r) and z max lie above the threshold for unusual origin. Note that the A component in 
NGC2451 is a nearby moving group and the B component a background object. 
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